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Abstract 



A general method is described for finding algebraic expressions for matrix elements of any one- and 
two-particle operator for an arbitrary number of subshells in an atomic configuration, requiring neither 
coefficients of fractional parentage nor unit tensors. It is based on the combination of second quantization 
in the coupled tensorial form, angular momentum theory in three spaces (orbital, spin and quasispin), and 
a generalized graphical technique. The latter allows us to calculate graphically the irreducible tensorial 
products of the second quantization operators and their commutators, and to formulate additional rules 
for operations with diagrams. The additional rules allow us to find graphically the normal form of the 
complicated tensorial products of the operators. All matrix elements (diagonal and non-diagonal with 
respect to configurations) differ only by the values of the projections of the quasispin momenta of separate 
shells and are expressed in terms of completely reduced matrix elements (in all three spaces) of the second 
quantization operators. As a result, it allows us to use standard quantities uniformly for both diagonal 
and off-diagonal matrix elements. 
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1 Introduction 



In order to obtain accurate values of atomic quantities it is necessary to account for relativistic and correla- 
tion effects. Relativistic effects may be taken into account as Breit-Pauli corrections or, in a fully relativistic 
approach, by starting with the Dirac-Coulomb Hamiltonian and wave functions defined in terms of four- 
component one-electron orbitals. In both cases, correlation effects may be considered either variationally 
or perturbatively. For complex atoms and ions, a considerable part of the effort must be devoted to cop- 
ing with integrations over spin-angular variables, occurring in the matrix elements of the operators under 
consideration. 

Many existing codes for integrating the spin-angular parts of matrix elements (Glass 1978, Glass and 
Hibbcrt 1978, Grant 1988, Burke et al 1994) are based on the computational scheme proposed by Fano (1965). 
In essence, it consists of evaluating recoupling matrices. Although such an approach uses Racah algebra, it 
may be necessary to carry out multiple summations over intermediate terms. Due to these summations and 
the complexity of the recoupling matrix itself, the associated computer codes become rather time consuming. 
A solution to this problem was found by Burke et al (1994). They tabulated separate standard parts of 
recoupling matrices along with coefficients of fractional parentage at the beginning of a calculation and used 
them further on to calculate the needed coefficients. Computer codes by Glass (1978), Glass and Hibbert 
(1978), Grant (1988), Burke et al (1994) utilize the program NJSYM (Burke 1970) or NJGRAF (Bar- 
Shalom and Klapisch 1988) for the calculation of recoupling matrices. Both are rather time consuming when 
calculating matrix elements of complex operators or electronic configurations with many open subshclls. 

In order to simplify the calculations, Cowan (1981) suggested that matrix elements be grouped into 
"Classes" (see Cowan 1981 Figure 13-5). Unfortunately, this approach was not generalized to all two- 
electron operators. Perhaps for this reason Cowan's approach is not widely used although the program itself, 
based on this approach is widely used. 

Many approaches for the calculation of spin-angular coefficients (Glass 1978, Glass and Hibbert 1978, 
Grant 1988, Burke et al 1994) are based on the usage of Racah algebra only on the level of coefficients 
of fractional parentage. A few authors (Jucys and Savukynas 1973, Cowan 1981) utilize the unit tensors, 
simplifying the calculations in this way, because use can be made of the tables of unit tensors and selection 
rules can be used to check whether the spin-angular coefficients are zero prior to computation. Moreover, 
the recoupling matrices themselves have a simpler form. Unfortunately, these ideas were applied only to 
diagonal matrix elements with respect to configurations, though Cowan (1981) suggested the usage of unit 
tensors for non-diagonal ones as well. 

All the above mentioned approaches were applied in the coordinate representation. The second quan- 
tization formalism (Judd 1967, Rudzikas and Kaniauskas 1984, Rudzikas 1991 and Rudzikas 1997) has a 
number of advantages compared to coordinate representation. First of all, it is much easier to find algebraic 
expressions for complex operators and their matrix elements, when relying on second quantization formal- 
ism. It has contributed significantly to the successful development of perturbation theory (see Lindgren and 
Morrison 1982, Merkclis et al 1985), and orthogonal operators (Uylings 1984), where three-particle operators 
already occur. Uylings (1992) suggested a fairly simple approach for dealing with separate cases of three- 
particle operators. Moreover, in the second quantization approach the quasispin formalism was efficiently 
developed by Rudzikas and Kaniauskas (1984). The main advantage of this approach is that applying the 
quasispin method for calculating the matrix elements of any operator, we can use the reduced coefficients of 
fractional parentage whose matrix elements are independent of the occupation number of the shell. All this 
enabled Merkelis and Gaigalas (1985) to work out a general perturbation theory approach for complex cases 
of several open shells. 
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Thus, it seems that it is possible to formulate an efficient and general approach for finding the spin- 
angular parts of matrix elements of atomic interactions, relying on the combination of the second quantization 
approach in the coupled tensorial form, the generalized graphical technique and angular momentum theory 
in orbital, spin and quasispin spaces as well as the symmetry properties of the quantities considered, which 
would be free of previous shortcomings. Gaigalas and Rudzikas (1996) suggested such an approach for finding 
matrix elements of any one- and two-particle atomic operator for the case of two open shells of equivalent 
electrons. But the situation is different when the matrix elements between more complex configurations are 
considered. 

An approach for the latter case is described in the present paper. One of the main ideas proposed here 
allows one to solve the problems related to the more complex configurations. Namely, we propose to apply 
Wick's theorem (see Lindgren and Morrison 1982) not in its usual general form while calculating the matrix 
elements, but rather only for groups of operators acting upon distinct shells of equivalent electrons. So, the 
ordering of operators generally would not be normal. 

It is a universal approach for finding algebraic expressions for matrix elements of any one- and two-particle 
operator in the general case of an arbitrary number of subshclls in an atomic configuration, heavily based on 
the exploitation of the quasispin technique and of the Wigner-Eckart theorem in quasispin space. Expressions 
for matrix elements and recoupling matrices are obtained by first classifying each matrix clement into one of 
four classes, depending on the number of subshells being acted upon. Each class is then further subdivided 
into cases and explicit expressions derived for each case in terms of triangular conditions, ry'-symbols, and 
reduced matrix elements. From these expressions, efficient procedures can be developed, that apply tests 
before performing computation. 

Tensorial expressions for any two-particle operator are presented in Section 2. They are based on the 
underlying assumption that the second quantization operators (both creation and annihilation), acting on 
the same open shell, must always be beside one another in a tensorial product and must be coupled into a 
resultant momentum. Then the second quantization operators, acting on the next shell, must follow, etc. 
Section 3 deals with the matrix elements between complex configurations. General expressions for recoupling 
matrices were found (Section 4) by use of the modified graphical technique of Jucys and Bandzaitis (1977), 
allowing us to calculate graphically the irreducible tensorial products of the second quantization operators 
and their commutators, and to formulate additional rules for operations with diagrams. The additional rules 
allow us to find graphically the normal form of the complicated tensorial products of the operators. All the 
graphical transformations we use here are fully described in Gaigalas and Rudzikas 1996. 

Exploitation of this new version of Racah algebra based on the angular momentum theory, on a generalized 
graphical approach, on quasispin approach, and on the use of reduced coefficients of fractional parentage 
for finding the spin-angular parts of two-particle operators is outlined in Sections 5-7. Some details of the 
calculations are presented in the appendix. 

2 Tensorial expressions for any two-particle operators 

In order to be able to find the expressions for matrix elements of the operators studied, we have to express 
these operators in terms of the irreducible tensors or their irreducible products. In this section we will 
present all the necessary tensorial expressions for any two-particle operator G. 

First, we express the operator in second-quantization form (Gaigalas and Rudzikas 1996) as 



G = 





(1) 



nili.rijlj.n'-l'n'.l'. 
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where, as is customary, the creation operators atdj appear to the left of the annihilation operators aj,aj, 
before defining the shells upon which the second quantization operators are acting. After defining the shells 
explicitly, the second quantization operators are transformed using their commutation relations so that all 
operators with the same nX (A = I, s) are beside one another. For example, in the case where the electron 
creation operator dj, and electron annihilation operators a], (where i = mlism^nis^ and a], act upon the 
same shell a, and operator aj acts upon another shell (3, we have: 



G (a[3aa) = \ ]T ]T a<t>± a™r'~l a% a '~) a&"± x 



x (n a A a mz a m Sa n / 3A / 3m^m Sfl |ff( KlK2fe ' CTlCT2fe Va^a^ a TO / Sa n a A Q mj' c( mg c( ) . 



(2) 



Here we imply that a tensorial structure indexed by {n\K 2 k, a\a 2 k) at g has rank «i for electron 1, rank k 2 
for electron 2, and a resulting rank k in the I space, and corresponding ranks o\<j 2 k in the s space. 

Now, applying the graphical approach of angular momentum theory (Gaigalas et al 1985), we can get an 
expression, namely 

G (af3aa) = \ ]T ]T ^_ 1 ^+a 12 - K ' 12 -a' 12 -k- P x 
x (n a X a npX p \\g^ k >°^\\n a X a n a X a ) x 

la la K 12 II s s °12 

K\ n 2 k > < (Ti u 2 k 

la 1/3 «12 J [ S S a 12 
1/2 ) lf3 la Ki 2 } \ S S CT 12 1 

12 k Ki \\ a' 12 k K s j 



1 /2 

X [Kl2,0-i 2 ] [k' 12 ,ct' 12 } 



> x 











(Kiifa)- 


a G/9«) x 


a ( '° s) x 




a x a 







(fcfc) 



(3) 



where [a, 6] = (2a + 1) (26 + 1), (noA^^A^! | 5 ( k i k 2*=,<ti<t2*:) | \n a X a n a X a ) is the two-electron submatrix (re- 
duced matrix) element of operator G and a is defined as (Judd 1967) 



/_^\/ a +s-m !ct -m s a t('a«) 



(4) 



nii a m s \ 1 —Tai ,—m s ' 

Expression (3) has summations over intermediate ranks k' 12 , a' 12 , Ki, K s in tensorial product. The angular 
momentum projection of (kk) is p, —p. 

In order to calculate the spin-angular part of a two-particle operator matrix element with an arbitrary 
number of open shells, it is necessary to consider all possible distributions of subshells, upon which the second 
quantization operators are acting. These are presented in Table 1. We point out that for distributions 2-5 
and 19 - 42 the shells' sequence numbers a, /?, 7, 5 (in bra and ket functions of a submatrix element) satisfy 
the condition a < (3 < 7 < S, while for distributions 6 - 18 no conditions upon a, (3, 7, 5 are imposed. 
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Table 1 : Distributions of subshells, upon which the second quantization operators are acting, that appear in 
the submatrix elements of any two-particle operator, when bra and kct functions have u open subshells 



No. 


di 


a 3 


4 


4 


submatrix element 


1. 


a 


a 


a 


a 


{...nj^ ...\\G{n t l l n l n%n'A l M...nJ^ ...) 


2. 
3. 
4. 
5. 


a 

a 

S3 


(3 
a 
(3 
a 


a 

S3 
S3 
a 


(3 
a 
a 

S3 


(...nJ^...n l N /...\\ 
^(nilinjl-jn'^n'jl'j) 
\\...nj^...n^...) 


6. 


a 


a 


S3 


S3 


(...n Q Z^...n^...||G||...n Q e^ 2 ---V? +2 ---) 


7. 
8. 
9. 
10. 


P 
a 

a 
a 


a 

13 
(3 
(3 


a 
a 

S3 
a 


a 
a 
a 

S3 


{...n a l^...n fj l^ 13 ...\\ 
Giniknjljn'^n 1 ^) 
\\...n a l^+\..n p lf-\..) 


11. 
12. 
13. 
14. 


(3 

i 
i 

(3 


7 
(3 
(3 
7 


a 

7 
a 

7 


1 
a 

1 
a 


{...nj^npl^n^'' ...\\ 
Giniknjljn'^n'/j) 


15. 
16. 


7 
7 


7 
7 


a 

S3 


(3 
a 


(...nj^npl^n^' 1 ... WG^iknjljn'^n'^) 
||...nX" +1 V^ + V^ T ~ 2 ---) 


17. 
18. 


a 

13 


(3 
a 


7 
7 


1 
1 


(...nal^npl^n^ly"' ... WG^iknjljn'^n'^) 
\\...n a l^- 1 n l^~ 1 n 1 l^'' +2 ...) 


19. 
20. 
21. 
22. 


a 

(3 
a 

13 


S3 
a 
(3 
a 


7 
7 
S 

s 


5 
5 

7 
7 


(nj^nplp 13 'n^'' 'n s lf s \\ 
^(nilinjljn'J^n'jl'j) 
\\nJ^- 1 n p l N /- 1 n 1 l^ +1 n s l^ +1 ) 
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Table 1 (continued) 



No. 






4 


4 


submatrix element 


23. 
24. 
25. 
26. 


7 
7 
<5 
<5 


S 
S 

7 
7 


a 

/? 
a 

/? 


p 

a 

P 
a 


Ginihnjljn&n'/j) 


27. 
28. 
29. 
30. 


a 
a 

7 
7 


7 
7 
a 
a 


/? 
/? 


S 

P 
P 
S 


{n a l^°'nf3l^' 3 rL 1 l^'' nslg s \\ 
G^iknjljn'^n'^) 

IKe-V^ + V^~Wf +1 ) 


31. 
32. 
33. 
34. 


/3 
(5 

P 
S 


p 

8 

P 


a 
7 
7 
a 


7 
a 

a 


(n a l^npl N ^n^n s l^\\ 
GimknjlXl'.n'/j) 


35. 
36. 
37. 
38. 


a 
S 
a 
5 


6 
a 
S 
a 


/? 
7 
7 
/? 


7 
P 
P 
7 


G^iknjljn'^n'jl'j) 


39. 
40. 
41. 
42. 


P 
7 

p 

7 


7 
P 
7 
P 


a 

S 
a 


S 
a 
a 
5 


{n a l^"npl^ l3 rL 1 l^''rL 8 l^ s \\ 
Giniknjljn&n'/j) 
||nX" + V^~V^~Wf +1 ) 



Let H be an array of intermediate coupling parameters in tensorial form, including K12, ovi, k' 12 , <j' 12 and 
possibly also others. Then the tensorial expressions for all these distributions can be grouped into four 
classes, where in each class, the two-particle operator G, operating on specific shells (see Eq. (1)), has one 
of four forms: 

1. All the second quantization operators act upon the same shell (distribution 1) and 



K12,0"12,k' 12 ,<tJ 2 P 

2. The second quantization operators act upon the two different shells (distributions 2-10) and 



(5) 



Kl2,Cl2,«^ 2 ' Cr i2 ^ 



(fcfc) 



p,-p 



(6) 



3. The second quantization operators act upon three shells (distributions 11-18) 
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G{III)~ ^6 (n Q A a ,n (3 A /3 ,n 7 A 7 ,S) x 

1 ( kk ) 



p,-p 



(7) 



4. The second quantization operators act upon four shells (distributions 19-42) and 



G(IV) 



K12,Ci2,k' 12 ,<tJ 2 P 
(K12CT12) 



(n Q A a , n^A^, n 7 A 7 , n^Aa, S) x 



(kk) 



(8) 



In (5)-(8), (nA, S),...,0 (n Q A a , npXp, n 7 A 7 , n^Aa, S) are proportional to the radial part of the operator 
G, and A^ kkS> (n\ 7 'E,),...,E( kk ) (nA, S) denote tensorial products of irreducible tensors. Parameter EE implies 
the array of coupling parameters that connect to the tensorial part. The explicit tensorial expressions are 
presented in the Appendix, using the graphical approach of Gaigalas et al (1985). Graphical methods make 
it possible to reduce the number of expressions from 42 to 6 for all distributions presented in Table 1 . Such a 
joining up of several distributions is possible by graphical means, because in the graphical technique of Jucys 
and Bandzaitis (1977), as in the tensorial products of operators of second quantization, the main elements 
are the Clebsch-Gordan coefficients. Therefore we may join up all the distributions having essentially the 
same algebraic structure, although with different tensorial products. The latter are represented by diagrams 
in which all the pecularities of a tensorial product are seen, and the differences of particular distributions 
are easily noticed. The use of other graphical methods (see e.g. Yutsis et al 1962 or Lindgrcn and Morrison 
1982) in joining up the distributions is complicated, since there the Wigner coefficients play the main role, 
and these are not fully compatible with the graphical transformations of the operators of second quantization 
in coupled form. Having classified the operators, we will now consider matrix elements of these operators 
for arbitrary configurations. 



3 Matrix elements between complex configurations 

Now, having the irreducible tensorial form of the operator being considered, we are in a position to find their 
matrix elements and recoupling matrices. Suppose that we have a bra function with u shells in LS coupling: 

V£™ (LSM L M S ) ee 

= (n 1 l 1 n2l2-n u l u a 1 L 1 S 1 Q 1 MQ 1 a2L 2 S2Q2MQ 2 ...a u L u S u Q u MQ u ALSM L Ms\ (9) 
and a ket function: 

^ {L'S'M' L M' S ) ee 

ee \n 1 l 1 n 2 l 2 ...nJ u a[L[S[Q[M^y 2 L 2 S 2 Q 2 M^ 2 ...a' u ^ (10) 

where A stands for all intermediate quantum numbers, depending on the order of coupling of momenta LiSi. 
Label Qi is the quasispin momentum of the shell riilf\ which is related to the seniority quantum number 
v>i, namely, Qi = (2Zj + 1 — Vj) /2, and its projection, MQ i = (Ni — 2k — 1) /2. In (9) and (10) on denotes 
all additional quantum numbers needed for the classification of the energy levels of the relevant shell. 

Using the Wigner-Eckart theorem in LS space we shift from the matrix element of any two-particle 
operator G between functions (9) and (10) to the submatrix element (V't™ (LS) ||G||^ e ' (L'S')) of this 
operator. 
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A general expression for the submatrix element of any two-particle operator between functions (9) and 
(10) with u open shells can be written as 



(^(LS)\\G\\^(L'S')) = 

E E E(- 1 ) Ae 'M-^' A ^ n ^'^ A ^ s ) x 

xT (mXi, n 3 X 37 n' l X' l7 n' J X' J ,A bra 7 A ket ,E, r) x 

xi?(A 4 ,A„A^.,A"™,A fcet ,r) (11) 

where A bra = {L i S i ,L j S j ,L' i S' i ,L' j S' j ) bra is the array for the bra function shells' terms, and similarly for 
A ket . This expression is similar to Eq. (136) used by Grant (1988) in his derivation. So, to calculate the 
spin-angular part of a submatrix clement, one has to compute: 

1. The recoupling matrix R (A i; X 3 ;, A^, A^, A bra , A ket , r) . This recoupling matrix accounts for the change 
in going from matrix element (^>„ (LS) \\G\\ip ket (L'S')) , which has u open shells in the bra and 
ket functions, to the submatrix element T {riiXi, n 3 X 3 , n\X^, n' 3 X^ A bra , A ket , S, T) , which has only the 
shells being acted upon by the two-particle operator in its bra and ket functions. 

2. The submatrix element T (n^Ai, n 3 ,X 3 , n-A-, n'jX'j, A 6ra , A ket , S, r) , which denotes the submatrix ele- 
ments of operators of the types A( kk ') (nA.S), B^ kk ') (nX, S), c( fcfc ') (nA, H), L>( Zs ), £( fefc ') (nA, S) (see 
(5)- (8)). Here V refers to the array of coupling parameters connecting the recoupling matrix 

R (A;, Xj, X'i, X' 3 ,A bra , A ket ,T) to the submatrix element. 

3. Phase factor A. 

4. 6' (niXi,n 3 X 3 ,n' i X l i ,n' 1 X' J ,'E)j, which is proportional to the radial part and corresponds to one of 
6 (nA, EE),.. .,6 (n a X a , npXp, n 7 A 7 , nsXs, 5). It consists of a submatrix element 

(niXin 3 X 3 \\g( KlK2k ' aia2k ^\\n' i X' i n J X J ) , and in some cases of simple factors and 3nj-coefficients. For 
instance, for the distributions aa/3/3, ~fja/3, 77/3a, a/377, P a llj fiaSj, a/3Sj, Pajd, jSafi, 

5 j (3a, 7<5/3a, 5ja(3 (see expressions (52), (53) and notes on 6' and in Appendix) it is: 

6' (n i Xi,n 3 X 3 ,n' i X[,n J X' j ,'EL) = (-1)* 6 (rijAj, n 3 Xj,n' i X' i ,n ] X J ,'EL) = 
= 1 (^AinjAjll^^^^^Hn^ni.Aj) x 



X [«12, CT12, «/ 12 , cr 12 ] 



1/2 



Kl2 «/ 12 




(12) 



where the integer £ determining the phase depends upon the configuration states involved. Rules for 
its determination are given in the Appendix. 

The calculation of 9' (njAj, n 3 X 3 , n^A-, n'jX' 3 ,E) is straightforward (from an angular momentum point of 
view) and depends on the radial form of the operator. In the next sections we will describe expressions for 
the recoupling matrix, the submatrix elements, and the phase factor, respectively. 



4 Recoupling matrices 

In this section we present the expressions for the recoupling matrices 

R (Aj, Xj, A-, X'j, A bra , A ket ,T) . 
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These matrices may be treated in the orbital I and spin s spaces separately. That is, 



R(X i ,X j ,\' i ,X' j ,A bra ,A ket ,T) = 
= R (l i ,l j ,l' i ,l' j ,A b l ™, Af et , r,) R (s, s, s, s, A bra , A^,T S ) 



(13) 



where A bra 



[Li, Lj,L[, L'j) bra and A b s ra = (Sj, Sj, Sj) bra . Therefore for simplicity we present only the 



expressions in I space. The recoupling matrices in s space are easily obtained from analogous expressions in 
I space by making corresponding substitutions l\, l 2 ,---,l u — * s; L\ — ► Si, L 2 — > S 2 ;...; ^12 — * 5*12, 
L123..U-1 — * S'i23...u-i; L — ► S, L' — ► S". Also, the analytical expressions for recoupling matrices 
presented in this section are valid in the case of jj-coupling. 

As we have mentioned earlier, there are four classes as defined by Eq.'s (5)- (8). We will consider each class 
separately. All the expressions presented below are obtained by using the approach of angular momentum 
theory described by Jucys and Bandzaitis (1977). 



4.1 One interacting shell 

Let us assume that the operators of second quantization act upon shell a as in distribution 1 of Table 1, 
where a = a. Then the recoupling matrix has the expression: 



R{l ai L ai k) 



[L a }- 1/2 S (L 1} Li) ...5 (L a _ 1; L;_0 S (La+i, i'„+i) -S (L u , L'J 



x < 



S(L u L[,k); 
Ci; 

C 1 C 2 (k,a + l,u-l)C 3 ; 

6 (L12, L' 12 ) ...S (Li2...a-1, i'l2...o-l) X 

xCiC 2 (fe,a+ 1)C 3 ; 



for m = 1 
for tt = 2 
for a < 3, u > 2 

for a>3, a ^ u,u > 2 



(14) 



5(Ll2,L' 12 ) — (-kl2...a-l,£'l2...a-l) ^3 ; fora = u,«>2 

In the above, the notation 5 (Li, L\, k) means the triangular condition | L\ — L[ |< k < L\ + L[ and 



a = [l ,t 



n i/2 \ k L' a L a 
I J T T 



(15) 



where the values of parameters <p, J, T and T" present in expression (15) are given in Table 2. The remaining 
two coefficients arc 



c 2 (k, fc min , fc max ) = n [Li2...i-i,Li 2 ...j] 

^rnin 



1/2 



k ^'l2...i-l ^12...»-1 



Li 



12.. 



and 



(16) 



(17) 
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Table 2: Parameters for equation (15) 



u 


a 




J 


T 


V 


2 


1 


Li + 2L[ -L 2 -L' + k 


Lz, 


L 


L' 


2 


2 


L x + L + L' 2 + k 


U 


L 


V 


u 2 


1 


L x + 2L[ -L 2 - L' 12 + k 


Lz, 


L\ 2 


n 12 


u ^ 2 


2 


Li + L 12 + L' 2 + k 


Li 


L\ 2 


T< 
n 12 


M ^ 2 


a > 2 


Ll2...a-1 + Li2... a + L' a + k 


L\ 2 ...a-l 


L\ 2 ... a 


L'l 2 ...a 



Table 3: Parameters for equation (17) 



u 




3 


J 


J' 


T 


T 


u ^ a 


fc + 


L u + L12...U-1 + L' 


L u 


Ll2... u -1 


^12...M-l 


L 


V 


a 


k- 


L\ 2 ... u -i + 2L U + L' u — L 


L\ 2 ... u -i 




L'u 


L 


L' 



11 



where the parameters ip, j, J, J' , T and T" are given in Table 3. 

When the total rank k = 0, the recoupling matrix becomes simply 

R (l a , L a , 0) = 6 (L U L[) S (L 2 , L' 2 ) S (L 12 ,L[ 2 ) ...5 (L a - U L' a _ x ) x 

XS (il2...o-l, i'l2...o-l) 5 ( L «> L 'a) 6 (£-12...a, L[ 2 J 5 (L a+1 , L' a+1 ) X 

xS (L 12 ... a+1 , L' 12 . a+1 ) ...6 (L u , L' u ) 8 (L, L>) . 
Expression (18) is equivalent to (13.60) of Cowan (1981). 



(18) 



4.2 Two interacting shells 

In this case let us assume that the operators of second quantization act upon the shells a and b (distributions 
2-10 in Table 1, where for distributions 2-5 a = a, b = (3 and for others (6-10) a = min{a, (3}, b = max{a, /?}). 
Then 

R (l a , L a , h, L b , K12, k' 12 , *) = (-1) C [La, Lb}^ 2 6 (L X , L[) ...S (L a - U L' a _^ X 

x<J (L a+1 ,L' a+1 ) ...6 (L b -i, L' b _ 1 ) 8 (L h+1 ,L' b+1 ) ...5(L U ,L' U ) x 



where 



and 



x < 



C 4 (K 12 ,K' 12 ,k,l)C 2 (k,3,u-l)C 3 ; 

C\C 2 (K 12 ,a + 1, b - 1) C 4 (K 12 , K[ 2 , k, 1) x 
xC2(fc,6 + l,«-l)C 3 ; 

C X C 2 (K 12 , a + 1, b - 1) C 4 (K 12 ,K' 12 , k, 1) ; 

S (L 12 , L' 12 ) ...S (L 12 ...a-l, ii 2 ... a _i) Ci x 
xC 2 (K 12 ,a + 1, b - 1) C 4 (K 12 , K[ 2 , k, 1) x 
xC2(fc,6 + l,«-l)C 3 ; 

S(L 12 ,L' 12 )...S(L 

12...a-l, L' 12 a _ 

i) Cix 

xC 2 (A-ia, a + 1, 6 - 1) C 4 (^12, JT( 2) k, 1) ; 



for a = 1,6 = 2 



for a < 3, 6 > 2, 6 ^ u 



for a < 3, 6 = u 



for a > 3, b > 2, b ^ u 



for a > 3, b - 



for a < /?, 

t!2 + «'i2 — for a > /?, 



C 4 (h,k 2 ,k,P) = [.h,J 2 ,J'z,k] 



1/2 



J( fcl Jl 
J 2 fc 2 J 2 

Jk k ,h 



> . 



(19) 



(20) 



(21) 



The values of parameters J\, J[, J 2 , J 2 , J3 and J3 present in the expression (21) must be taken from 
Table 4. For the case a < (3 in Eq. (19) K 12 = m 2 , K' 12 = k' 12 and when a > (3, then K 12 — k' 12 , K[ 2 = m 2 . 
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Tabic 4: Parameters for equation (21) 



p 


a 


b 


u 


Ji 


J'l 


J 2 


4 


J 3 


4 


1 


1 


2 


u ± b 


Li 


L[ 


i 2 


L' 2 


L\2 




1 


1 


2 


b 


Li 


L[ 


i 2 


L' 2 


L 


L 


1 


a 


+ \ b^2 


b 


L\... u -\ 




L u 


L' u 


L 


V 


1 




in all other 


cases 


L\...b-\ 


^l...b-l 


L b 


L' b 


Ll...b 


L[... b 


2 


in all cases 


L\... c -i 


ii...c-l 


L c 




L\... c 





When the total rank k = 0, and k\ 2 = k' 12 — the recoupling matrix has the form: 

R (l a ,L a , l b , L b , k, k, 0) - [L a , L' b , k}- 1/2 S {L u L[) ...S (L a _ u x 
x<* (L a+1 , L' a+1 ) ...5 (Lb-uL^) S (L b+1 ,L' b+1 ) ...S {L u , L' u ) x 
xS (L 12 ,L[ 2 ) ...S (L 12 ... a -u i'l2...a-l) s {Li 2 ... b , L' 12 _ b ) ...S (L, L') X 

C 5 (1) ; for a = 1,6= 2 

CiC 2 (M+l,&-l)C 6 (l); 

X < 

5 (Ll2, i' 12 ) -5 (L 12 ... a -U L' 12 ...a-l) X 

k xCiC 2 (fc,a+l,6-l)C 5 (l); 

where 

C 5 (P) = (-l) k+L » +J ^[JuL>} 1/2 { * % ^j. (23) 

The values of parameters Ji, J{ and J 2 present in the expression (23) must be taken from Table 5. 

Formula (22) has no analogue in Cowan (1981). Our expressions for the recoupling matrix do not depend 
on coefficients of fractional parentage and have no intermediate summations. Therefore they will be very 
convenient for practical calculations. 

4.3 Three interacting shells 

When the operators of second quantization act upon three shells a, b and c (distributions 11-18 in the Table 
1), we have: 



for a < 3 
for a > 3 
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Tabic 5: Parameters for equation (23) 



p 


Case 


Ji 


A 


J2 


1 
1 
1 


a = 1 and 6 = 2 
b = u 


L a 

Li...b-i 
Li...b-i 


L' a 
Jl 

^1...6-1 


L\2...b 
L\2...b 

L 


2 
2 


c^u 
c~u 


L\... c -\ 
L\... c -i 


Jl 

M...C-1 


L\2...c 

L 


3 
3 


d^u 
d = u 


L\...d-i 
L\...d-i 


Jl 
Jl 

L \...d-\ 


L\2...d 

L 



R (la, L a , lb, Lb, l c , L c , k\, k2, Kl2, k'i 2 , — 

= [L a , L b , L c ]- 1/2 5 (L 1} Li) ...5 (L a _ 1; L^) x 
:S (L a+l ,L' a+1 ) ...8 (L 6 _i, L^) 5 (L 6+1 , L' 6+1 ) ...5 (L c _ 1; L^) x 
x<5 (L C+1 ,L' C+1 ) ...5 (L u , L' u ) J2 (-1) C C 6 x 

312 



X < 



C 4 [K u K2,hi, 1) C 2 Cm, 3, c - 1) x 
x C 4 (ji2 , if 3, k, 2) C 2 (A, c + 1 , u - 1) C 3 ; 

CiC 2 (ifi, a + 1, 6 - 1) d (if ! , K 2 ,j 12 , 1) x 
XC 2 (jl2, b + 1, c - 1) C 4 (ji2, if 3, fc, 2) X 
xC 2 (fc,c+l,«-l)C 3 ; 



for a = 1,6= 2 



for a < 3 



for a > 3, 



<5(Li 2 ,L' 12 )...S(L 

12...a-l, L' 12 ... a _i )Cix 

xC 2 (if l5 a + 1, 6 - 1) C 4 (if i, if 2 ,. 7i2, 1) x 
xC 2 (ii 2 , 6 + 1, c - 1) C 4 (ji 2 , if 3 , fc, 2) x 
xC 2 (k,c+l,u-l)C 3 ; 

where parameters a, b, c, C, ifi, if 2, if 3 an d coefficient Cq are given in Table 6. The coefficient 
C' 6 (ki,k 2 ,k 3 ,k4,k 5 ,k e ) is 



ki k 2 k; 
fc 4 

From (7) we have that in expressions (24) and (25) the ranks fc 4 = l a , k 2 = lp. 



Cafcl,fc 2 ,fc3^4,fc5,fc 6 ) = (-l) fel + fe2 - fe3+2fc5 [fc3,fc 6 ] 1/2 <| I' I' I" 

k 4 k 5 k e 
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Tabic 6: Parameters for equation (24) 



Case 


a 


b 


c 


c 


Ki 


K 2 


K 3 


c 6 


a < (3 < 7 


a 


P 


1 





ki 


k 2 


K 12 


5 (jl2, K12) 


P < a < 7 


P 


a 


1 


kl + k>2 - Kl2 


k 2 


ki 


K 12 


<5 (jl2, K12) 


P < 7 < a 


P 


1 


a 





k 2 


K 12 


ki 


C 6 («'l2>*2,jl2,fel,A;, K12) 


a < 7 < P 


a 


1 


P 


ki + k 2 - K12 


ki 


K 12 


k 2 


C<3 Wi 2 ,ki,ji2,k 2 ,k,Ki 2 ) 


7 < a < P 


1 


a 


P 


2ki +k2- K12+ 


K 12 


ki 


k 2 


Cg (k'i2,A;i, ji2,k 2 ,k, K12) 










+«'l2 - 312 










7 < P < a 


1 


P 


a 


ki + k 2 - K12 


K 12 


k 2 


ki 


C<3 («'l2>* ; 2,jl2,fcl,A;,Kl2) 



When the total rank fc = 0, and K12 = k/ 12 = fc, the recoupling matrix has the form: 

R(l a ,L a ,l b ,L b ,l c ,L c ,ki,k 2 ,k,k,0) = (-l) c [L a , £ 6 , L^,, i^p 172 x 
x<f (£1, Li) ...5 5 (L„+i, L; +1 ) ...<5 (L 6 _x, L'^) x 

xJ (L 6 +i, ^6+1) ...<5 (Lc-i, ic-i) * (^+1-^+1) •••< 5 L «) x 

X<5 (L12, L' 12 ) ...<5 (Ll2...a-1, L'l2...a-l) 5 (^12...c, A 2 ...c) ■•■<* (£> £') X 



d (Ki,K 2 ,K 3 , 1) C 2 (^ 3 , b + 1, c - 1) C 5 (2) ; for a = 1, & = 2 



x < 



dC 2 (Ki, a + 1, b - 1) C 4 (#i, # 2 , #3, 1) x 
xC 2 (K3,&+l,c-l)C 5 (2); 

(5 (L12, L' 12 ) (Ll2...a-1, L'l2...a-l) Cl X 

xC 2 {K u a + 1, 6 - 1) C 4 (Ki,K 2 , #3, 1) x 
[ xC 2 (]f 3 ,6+l, C -l)C 5 (2); 



for a < 3 



for a > 3 



(26) 



where the parameters £, Ki, K 2 , K 3 values are given in Table 7. 

The recoupling matrix for three interacting shells (26) has the same advantages as the equivalent quantity, 
Eq. (22), for two shells. 



4.4 Four interacting shells 

When the operators of second quantization act upon four shells a, b, c and d (distributions 19-42 in the 
Table 1), we have: 
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Tabic 7: Parameters for equation (26) 



Case 


c 




K 2 


^3 


a < (3 < 7 





ki 


k 2 


fc 


P < a < 7 


fci + k 2 - fc 


k 2 


fci 


fc 


P < 7 < a 


2fci 


k 2 


k 


fcl 


a < 7 < p 


k\—k 2 — k 


ki 


k 


fc 2 


7 < a < P 


2k 


k 


fci 


fc 2 


7 < P < a 


ki + k 2 + k 


k 


fc 2 


fcl 



i? (Z a , L a , h, Lb, l c , L c , Id, Ld, fci, fc 2 , Ki 2 ,k 3 ,k 4 , k' 12 , fc) = 
= [L a , L b , L c , L d p 1/2 5 (Li, Li) ...(5 (L a _i, L^) x 
?5 (L a+1 , L' a+1 ) ...(5 (L b _i, L b _ 1 ) <5 (L 6+1 , L' b+1 ) ...J (L c _i, L^) x 
xS (L C+1 ,L' C+1 ) ...8 (L d _i, L' d _^) (5 (L d+1 , L^ +1 ) ...(5 (L„, L^) x 



C 4 (fcl,fc2,Kl2,l)C2(Kl2,3, c- l)C 7 (c,rf) x 

xC 2 (fc,rf+l,u- 1)C 3 ; 

CiC 2 (fci,a+ 1,6-1) C 4 (fci,fc 2 , «i2, 1) x 
xC 2 (ni 2 ,b+ l,c- 1) C 7 (c, d) x 
xC 2 (fc,rf+l,u-l)C 3 ; 

S(L 12 ,L[ 2 ) ...6 (L 12 ... _i,Li 2 .. _ 1 ) Cix 
xC 2 (fci, a + 1, b - 1) C 4 (fci, fc 2 , K12, 1) x 
xC 2 (ki 2 ,6+ l,c- 1) C 7 (c, <f) x 
xC 2 (fc,rf+l,u-l)C 3 ; 



for a = 1,6 = 2 



for a < 3 



for a > 3, 



C7 (fcmin, fcmax) 

^C 8 (I) C10 (/) , for fc max - fcmin = 1 

max ^min 2 



2EC8(Ii)Cii(Ji,/2)Cio (J 2 ); 



for /c max &min ^ 2 
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Cs(I) = (-l)"»+ i '»-c- J [L C) /,Li2...c-l,Li2...c] 



fc 3 

Ll2... c -1 



1) K12 


+i'l2... 




L c 


Ll2.. 





1/2 



12...C-1 K 12 
Ll 



C 9 (/i,/ 2 ,i) = (_i) 2 ^+ i *)+ii 2 ...i+ii 2 ... 1 +fcs+i«i a [£i2 ___._ 1}Jl) 7 2 ,L' 12 .J 



1/2 



Ll2...»-1 ^1 
Lj Ll2.. 



fc 3 



7', 



12.. .i-l 
/ 2 



L' 



12. ..i 



«12 
Li 



(29) 



(30) 



ClO (-0 = ( _ l) 2 ( / + fe3 ) + ' i:4+Kl2 + K 'l2 + fe + i l 2 -- d + L 'l2...d+- L £i + i d+ i 'l2 ...d-1 x 

x [/sia, «i 2 , L d , I, L' 12d , ii 2 ... d -i] 1/2 ^ (-1)" [x] X 



I 7 k' 12 x 

\ k L 'l2...d-l «12 

L\2...d-\ h x 



^12 *^ 
&4 Ll2...d-1 h 



7', 



Cn (J 1 ,/ 2 ) = (-l) 



Ll2...d L d 
Ji-/ 2 +L' 15 



v 12...d-l 
Ll2...d 



T 1 T' 

L d h l2...d 

,1/2' 



xC 2 (i,c+l,d-l) 



-M/i,/ 2 ] 1/2 ENx 

f fc 3 



J 2 



fc 3 Kl2 X I I fc 3 Kl2 

v, L 12...c L 12...c h J } L'i2 ...d-1 L 12...d-1 

From (8) we have that in the expressions (27), (30), (31) and (32) the ranks k\ = l a 
fc 4 = h. 

When the total rank k = and «i 2 = k' 12 — k, the recoupling matrix has the form: 

R (l a , L a , l b , Lb, l c , L c , l d , L d , ki,k 2 ,k, k 3 , fc 4 , k, 0) = 
= [L a , L h , L c , L' d , k}- 1/2 6 (Li, Li) ...S (7 a _ 1; 7^) x 
X(5 (L a+1 ,L' a+1 ) ...5 {Lb-i, L' b _-y) 5 (L b+1 ,L' b+1 ) ...x 
x...S (Lc-uL'^j) 5 (L C+1 ,L' C+1 ) ...5 {L d -x,L' d _^) x 
x5(L d+1 ,L' d+1 ) ...S(L U ,L' U ) x 
xS(L 12 ,L[ 2 )...S(L 

12...o— 1) L 12 a _ l )6(L 12 ...d,L> 12 ^ d )...6(L,L>)x 



1/3, h 



(31) 

(32) 
— I j, 



C i (k 1 ,k 2 ,k,l)C 2 (k,b+l,c-l)x 
xC 4 (fc,fc 3 ,fc 4 ,l)C 2 (fc 4 ,c+l,d- 1) x C 5 (3) ; 

CiC 2 (fc, a + 1, 6 - 1) C 4 (fei, fc 2 , fc, 1) x 
xC 2 (fc, 6 + 1, c - 1) C 4 (fc, fc 3 , fc 4 , 2) x 
xC 2 (fc 4 ,c+l,rf- 1)C 5 (3); 



for a = 1,6 = 2 



for a < 3 



(33) 



for a > 3 



(5 (L12, L' 12 ) ...<5 (Li2... a _i, L' 12 . „ a _i) Ci x 
xC 2 (fc, a + 1, 6 - 1) C 4 (fci, fc 2 , fc, 1) x 
xC" 2 (fc, 6 + 1, c - 1) C 4 (fc, fc 3 , fc 4 , 2) x 
xC 2 (fc 4 ,c+M-l)C 5 (3). 

Expression (33) have also no analogue in Cowan (1981). 

Thus, we have studied all possible cases of matrix elements of arbitrary two-electron operators. The 
expressions for recoupling matrices ((22), (26) and (33)) obtained in this section are simpler and, thus, more 
convenient for practical applications, than those of Cowan (1981), except for the simplest case fc = of an 
operator acting on one shell (18), where they are equivalent. 
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5 Calculation of tensorial quantities 



In this section we will consider the submatrix elements 

T (riiXi, nj Xj ,n'X, »J X' j ,A bra , A ket , S, r) 



appearing in (11). Taking into account, that operators am\ an d a m ' x are components of the tensor cim 



,(9A) 



q m x , 



having in quasispin space the rank q = | and projections m q = ±5, i.e. = am}m s and a^i 



,(?A) 



the operators A^ kk ') (nA.H), s( fefc ') (nA, S), c( fcfc ') (nA, S), £><'">, ^( fefc ') (nA, S) (sec (5)-(8)) in our case cor- 
respond, respectively, to the following five expressions: 



j (9 A) 



(kicti) 



(kicti) 



(«i<7i) 



x a 



(gA) 

"1,3 



(K2CT2) 



(K2C2) 



,(?A) 



,(«A) 
"m g2 



(kicti) 



,(?A) 
m 9 3 



,(gA) 

-"771,4 



(K2CT2) 



-1 (fefc) 



(34) 
(35) 

(36) 

(37) 

(38) 



We will discuss the derivation of submatrix elements of these operators, and present the expressions for 
these quantities. It is worth noting that these tensorial quantities all act upon the same shell. So, all the 
advantages of tensor algebra and the quasispin formalism may be exploited efficiently. 

We obtain the submatrix elements of operator (34) by straightforwardly using the Wigner-Eckart theorem 
in quasispin space: 



(l N a QLS\\a^: ] \\l N ' a'Q'L'S 1 ) = - \Q]~ 



1/2 



x (I aQLS\\\a^\\\la'Q'L'S') 



Q' 1/2 
M' Q m q 



M Q 



(39) 



where the last multiplier in (39) is the so-called completely reduced (reduced in the quasispin, orbital and 
spin spaces) matrix element. The coefficient 



3i 32 3 
mi TO2 m 



is a Clebsch-Gordan coefficient. Different notations for it occur, e.g. A^^ m2 in Eckart (1930), Sjl^ 2 im2 in 
Wigner (1931), (jij2'mi'm,2\jij2jm) in Condon and Shortley (1935) or Judd (1967). 

The value of the submatrix element of operator (35) is obtained by basing our development on (33), (34) 
of Gaigalas and Rudzikas (1996). In the other three cases (36), (37), (38) we obtain them by using (2.28) of 



18 



Jucys and Savukynas (1973): 

{nl N aQLS\\ [F^l (nX) x (nX)] (kk) \\nl N ' a'Q'L'S') = 

x {nl N aQLS\\F {K ^ {n\)\\nl N " a"Q"L"S")x 

a"Q"L"S" 

x(nl N " a"Q"L"S"\\G^^ (nX) \\nl N ' a'Q'L'S')x 

X \L< L L»)[s> S (40) 

is one of (34) or (35) and the submatrix elements correspondingly are defined 
by (39) and (33), (34) of Gaigalas and Rudzikas (1996). N" is defined by the second quantization operators 
occurring in F (Kiai '> (nX) and G (K2CT2) (nX). 

As is seen, by using this approach, the calculation of the angular parts of matrix elements between 
functions with u open shells is reduced to requiring the submatrix elements of tensors (34) and (35) within 
one shell of equivalent electrons. As these completely reduced submatrix elements do not depend on the 
occupation number of the shell, the tables for these quantities are considerably reduced in size in comparison 
with the tables of analogous submatrix elements of tensorial quantities U k , V klk2 (Jucys and Savukynas 
1973) and the tables of fractional parentage coefficients. 



6 Phase factor 

In this section we present the phase factors A in (11), which appear for submatrix elements of operators in 
Eq. (5)-(8). 

For distributions 1-6 (Table 1): 

A = 0. (41) 

For the distributions 7-18 (Table 1): 

i-i 

A = 1+ £ N k , (42) 

fe=» 

where if a < (3, then i = a, j = (3, and if a > {3, then i — (3, j = a; Nk is the occupation number of a shell 
of equivalent electrons having the label k. For the distributions 19-42 (Table 1): 

0-1 ,5-1 

A =]T N k + N k . (43) 

k— a k—~y 

7 Spin-angular part of any two-particle operator 
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Tabic 8: Scheme of the expressions for matrix elements of any two-particle operator 







G(T) 


















No. 


G 


a 





! 


s 


e 


R 




A 


j_. 


(AA\ 
l 44 J 












(A K\ 
(40) 


l i4 i, 




C41 1 




('47"! 




(6Q) 








(48), (4V) 


l 14 J, 


l i8 J 


(41 1 


2. 


(50) 


(6) 


(35) 


(35) 






(51) 


(19), 


(22) 


(41) 


3. 


(50) 


(6) 


(35) 


(35) 






(51) 


(19), 


^22) 


(41) 


4. 


(54) 


(o) 


(35) 


(35) 






(55) 


(Vd J, 


ZZ J 


(41) 


5. 


f54l 




(60) 


('35 N I 

(60) 






(00) 




'22) 


(411 


o . 


(UZ) 




(60) 


(60) 






(06) 


ri q) 

iiy ) , 


ZZ 1 


(41 1 


( . 


(00) 


(°) 


(60) 


(64) 






(06) 


C1 OA 


f99l 


(/191 
[41) 


fi 
o . 


(00) 




(60) 


(64) 






(OV) 


l ly 7, 


(991 




9. 


(60) 


(6) 


(34) 


(37) 






(62) 


(19), 


(22) 


(42) 


10. 


(60) 


(6) 


(34) 


(37) 






(63) 


(19), 


(22) 


(42) 


11. 


(50) 


(7) 


(34) 


(34) 


(35) 




(51) 


(24), 


(26) 


(42) 


12. 


(50) 


(7) 


(34) 


(34) 


(35) 




(51) 


(24), 


(26) 


(42) 


13. 


(54) 


(7) 


(34) 


(34) 


(35) 




(55) 


(24), 


(26) 


(42) 


14. 


(54) 


(7) 


(34) 


(34) 


(35) 




(55) 


(24), 


(26) 


(42) 


15. 


(52) 


(7) 


(34) 


(34) 


(35) 




(53) 


(24), 


(26) 


(42) 


16. 


(52) 


(7) 


(34) 


(34) 


(35) 




(53) 


(24), 


(26) 


(42) 


17. 


(52) 


(7) 


(34) 


(34) 


(35) 




(53) 


(24), 


(26) 


(42) 


18. 


(52) 


(7) 


(34) 


(34) 


(35) 




(53) 


(24), 


(26) 


(42) 


19. 


(52) 


(8) 


(34) 


(34) 


(34) 


(34) 


(53) 


(27), 


(33) 


(43) 


20. 


(52) 


(8) 


(34) 


(34) 


(34) 


(34) 


(53) 


(27), 


(33) 


(43) 


21. 


(52) 


(8) 


(34) 


(34) 


(34) 


(34) 


(53) 


(27), 


(33) 


(43) 


22. 


(52) 


(8) 


(34) 


(34) 


(34) 


(34) 


(53) 


(27), 


(33) 


(43) 
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Table 8 (continued) 







G(T) 


















No. 


G 


a 


(3 


7 


S 


e 




R 


A 


23. 


(52) 


(8) 


(34) 


(34) 


(34) 


(34) 


(53) 


(27) 


, (33) 


(43) 


24. 


(52) 


(8) 


(34) 


(34) 


(34) 


(34) 


(53) 


(27) 


, (33) 


(43) 


25. 


(52) 


(8) 


(34) 


(34) 


(34) 


(34) 


(53) 


(27) 


, (33) 


(43) 


26. 


(52) 


(8) 


(34) 


(34) 


(34) 


(34) 


(53) 


(27) 


, (33) 


(43) 


27. 


(50) 


(8) 


(34) 


(34) 


(34) 


(34) 


(51) 


(27) 


, (33) 


(43) 


28. 


(54) 


(8) 


(34) 


(34) 


(34) 


(34) 


(55) 


(27) 


, (33) 


(43) 


29. 


(50) 


(8) 


(34) 


(34) 


(34) 


(34) 


(51) 


(27) 


, (33) 
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In the previous sections, all the expressions required to calculate the spin-angular part of any two-particle 
operator were given. For convenience, the structure of the expressions (the numbers of the corresponding 
formulas) are summarized in the Table 8 for each distribution given in the Table 1. The classification numbers 
of the distributions are presented in the first column of the Table 8. The equation number of the tensorial 
expression of the two-particle operator G is given in the second column, and the equation number of the 
tensorial class of the two-particle operator, denoted by G(T), in the third column. 

In the next four columns are given the numbers of the formulas of the tensors, which act inside the shell. 
A tensor acting upon the a shell is given in the a column, and in the columns (3, 7, S - upon the /?, 7, S shells, 
respectively. Consequently, if we want to find submatrix element, T (n i X i7 njXj,n' i X' i ,njX'j 7 A bra ,A ket ,E,T), 
at first we have to calculate the submatrix clement of the tensor from column a between functions, consisting 
only of the a shell, and then to look for the submatrix element of the tensor from column (3 between functions, 
consisting only of /3 shell, and so on. Thus, we need to calculate only submatrix elements of the tensors 
acting upon a certain shell. The details of the calculation of these submatrix elements were discussed in the 
Section 5. 

The coefficients G are given in the column. The numbers of expressions for the recoupling matrix 
R (Aj, Xj, X' i7 X' 3 A bra , A ket , r) and phase factor A are given in the last two columns. From this table it is 
easy to derive the general formulas for spin-angular parts of matrix elements of any two-particle operator. 
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8 Conclusions 



The approach to matrix element evaluation that we present, is based on the the combination of the angular 
momentum theory as described in Jucys and Bandzaitis (1977), on the concept of irreducible tensorial sets 
(Judd 1967, Rudzikas and Kaniauskas 1984), on a generalized graphical approach (Gaigalas et al 1985), on 
the quasispin approach (Rudzikas and Kaniauskas 1984), and on the use of reduced coefficients of fractional 
parentage (Rudzikas 1991, Rudzikas 1997, Judd 1996). All this, in its entity, introduces a number of new 
features, in comparison with traditional approaches: 

(1) The tensorial expressions of a two-particle operator, presented in Section 2, allow one to exploit all 
the advantages of a new version of Racah algebra based on quasispin formalism when the latter is applied 
within each particular shell only. In particular, this is not only a reformulation of spin-angular calculations 
in terms of standard quantities, but also the determination beforehand from symmetry properties, of which 
matrix elements are equal to zero without performing further explicit calculations. That is determined from 
the submatrix elements T (rijAj, rijX } -,, n-A-, n'jX'j, A bra , A ket , S, V) . 

(2) It enables one to use the Wigner-Eckart theorem in quasispin space. This provides an opportunity 
to use tables of reduced coefficients of fractional parentage and tables of other standard quantities (Section 
5), which do not depend on the occupation number of a shell of equivalent electrons. Thus, the volume of 
tables of standard quantities is reduced considerably in comparison with the analogous tables of submatrix 
elements of tensorial operators U k , V kl and the tables of fractional parentage coefficients. This undoubtedly 
makes the inclusion of shells of equivalent / electrons with arbitrary occupation numbers considerably easier, 
and the process of selecting the standard quantities from the tables becomes simpler. 

(3) The tensorial form of any operator presented in Section 2 allows one to obtain simple expressions for 
the recoupling matrices (Section 4). Hence, the computer code based on this approach would use immedi- 
ately the analytical formulas for recoupling matrices R (Aj, Xj, X[, X' 3 ; , A bra , A ket , r) . This feature also saves 
computing time, because i) complex calculations leading finally to simple analytical expressions (Bar-Shalom 
and Klapisch 1988) are avoided, and ii) a number of momenta triads (triangular conditions) can be checked 
before the explicit calculation of a recoupling matrix leading to a zero value. These triangular conditions 
may be determined not only for the terms of shells that the operators of second quantization act upon, as 
is the case for the submatrix elements T (njAj, rijXj, n^A-, n'jX'j, A bra , A ket ,E, T) (see conclusion 1), but also 
for the rest of the shells and resulting terms. 

In this approach both diagonal and non-diagonal with respect to configurations, matrix elements are 
considered in a uniform way, and are expressed in terms of the same quantities. The difference is only in the 
values of the projections of the quasispin momenta of separate shells. 

In this paper all the expressions needed in the spin-angular parts of matrix elements of two-particle 
operators calculation are presented. This approach is applicable to one-particle operators as well. While 
calculating the spin-angular parts of the latter, all the expressions needed are included in the cases discussed 
for the two-particle operator. For instance, in the recoupling matrix calculation two of the four cases discussed 
above appear, namely, when all the second quantization operators act upon the same shell (Section 4.1) and 
when they act upon two different shells of equivalent electrons (Section 4.2). Thus, this approach is applicable 
to any one- and two-particle operator. Practical usage shows that a scries of difficulties persisting in the 
traditional approach to the calculation of angular parts of matrix elements based on the usage of coefficients of 
fractional parentage and unit tensors can be avoided and high efficiency may be achieved. Indeed, preliminary 
calculations show that computer programs based on our approach on average are 4-6 times faster than the 
other well-known codes (Gaigalas et al 1995). This methodology can easily be generalized to cover the case 
of relativistic operators and relativistic wave functions. 



22 



Acknowledgements 



This work is part of a co-operative research project funded by National Science Foundation under grant No. 
PHY-9501830 and by EURONET PECAM associated contract ERBCIPDCT 940025. One of the authors 
(CFF) was supported by a grant from the Division of Chemical Sciences, Office of Basic Energy Sciences, 
Office of Energy Research, U.S. Department of Energy. 



Appendix 

Here algebraic expressions are presented for the two-particle operator (1) in the irreducible tensorial form 
for all the distributions from Table 8. Although there are quite a few distributions, the structure of their 
algebraic formulas is similar, and therefore on the basis of a graphical approach (Gaigalas et al 1985) the 
expressions may be written down in a compact form, where one general formula includes all the cases having 
the same structure. Each particular formula is obtained from these by performing elementary graphical 
transformations according to the rules explained below. The general expressions are: 
(1) Distribution aaaa (case 1 from Table 8): 

For this distribution the analytical expressions (7), (8) in Gaigalas and Rudzikas (1996) are used, in 
which the quantum numbers njlj, n^, n'jl'j acquire particular values of a shell a. For the first form 
(Figure 1, A\), we have 



A 1 



Kl2Cl2^i 2 fT 12 



0/ (n a X a , n a X a , n a X a , n a X a , c) x 

-i (fefc) 



[ a (i-») x a <^)] (Kl2CTl2) x [a^ x fiC-*)]^' 12 "") 



(44) 



where 



0/ { n a^aj n a^cn n a^aj "a^ai ") — O ( n a^ai n a^ai n a^aj ''a^aj — 

= O (fl a X a ,E) 



(45) 



and 



0/ {jl a X ai Ti a X a: 7l a X a: 7l a X a: l^j) 
2 



1 ( -i xfc-p+1 r / / il/2 w 

= 2(-l) y [Ki2,(Ti2,K' 12 ,cr' 12 J X 



x {n a X a n a X a \\g (K ^ k ^ rj ^ 

( i i «•! i r 

t, a i a nj\ 
X < l a l a K 2 > < 
Kl2 K12 fc , 

In an equivalent second form (Figure 1, A 2 + A 3 ), we have 



s s o\ 

S S <7 2 



k 



(46) 



A 2 + A 3 = 0// a (n a X a ,n a X a ,n a X a ,n a X a ,E) x 
\ a (i a s) x ~ a (i a s)Y^) x r 5 (ia.) xa do 



i (K2CT2) 



(fcfc) 



p,-p , 
(fefe) 



+6//6 (n a A a ,n a A a ,n a A a ,n a A a ,S) [a('° s ) x a (iaS) ] p _ 



(47) 



where 



©//a (n a A Q ,r! Q A Q ,n Q A a ,n a A Q ,s) = \ (-l) fc p (Ji, k 2 , c 2 ] 1/2 x 



-1/2 



x in, 



A a n a X a \\g^ k ^^\\n a X a n a X a ) 



(48) 
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n n Xn 
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KK 



kk 
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riiXi 



n'X' 



kk 



rijXj riiXi 



— rr~- 

kk 



rijXj 
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K12CT12 



n'X' 



O 



1 

kk 



*12"12 



33 



riiXi 



33 



K12C12 K,,(7 



IT 

kk 



Y1 U Y1 



TljXj 



A 7 



As 



Figure 1: Diagrams for an arbitrary two-particle operator. Diagrams A\, A 2 and A 3 represent two-particle 
operators when this operator has distribution aaaa. Diagrams A 4 and A 5 represent two-particle operators 
for all others distributions. Diagrams A±, A 2 , A3, A 4 and A 5 arc similar to the usual Feynman-Goldstone 
diagrams. Diagrams Aq, A-j and As represent tensorial products of second quantization operators (A$ for 
the second group, Aj for the third, and A s for the fourth.). 
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and 



9/76 ( n a\i n a^aj "a^aj n a^a, 



-1) 



fe-p+1 



* A Q n Q A Q | l^i^^fc) | |„ Q A Q n Q X a ) x 



Kl K 2 fc 
l a l a l a 



g i a 2 k 

s s s 



(49) 



Both factors 6// a (n a X a , n a X a , n a X a , n a X a , S) and 0//;, (n a A Q , n a X a , n a X a , n a X a , S) have properties anal- 
ogous to those of 6/ as stated in (45). 

(2) Distributions a(3a(3, /3a(3a, P^a-f, j(3^a, a-fPS, 7a<5/3, (3Sa^ 7 S(3ja, aS/3j, Sajfi, (3^aS, r )f35a (cases 
2, 3, 11, 12, 27, 29, 31, 32, 35, 36, 39, 40 from Tabic 8) (Figure 1, A A , A 6 ): 



where 



A 4 = (niAi.rij-Aj-, n^,n^- A^-,S) A 6 , 

8 (nA^-A^n^n^-,-) = | (-l) fe ~ p [ki, cti, k 2 , er 2 p 1/2 x 
x (nAnAll3 (KlK2fe ' fflCT2fe) IK A H A j) ' 



(50) 



(51) 



Diagram A 6 corresponds to tensorial products of the operators of second quantization for two-particle oper- 
ator (for details see in Gaigalas and Rudzikas 1996). 

(3) Distributions aa/3/3, 770;/?, 77/3a, a/377, P a llj afi^fS, [3a5j, a/35j, PajS, 7<5a/3, 5jf3a, j5/3a, Sjaf3 
(cases 6, 15-26 from Tabic 8) (Figure 1, A 5 , A 7 ): 



A 5 = ^2 (niXi,njXj,n' i X' i ,n' j X' j ,'B) A 7 , 

^12^12^i 2 Cr 12 



(52) 



where 



fc-p+i 



6 (n i )^,n j X j ,n' i X' i ,n' j X' j ,'Z) = \ (-1) 
x[Ki 2 ,a 12 ,K' 12 ,a; 2 ] 1/2 (n J A J n J A J || 5 ^^ fc ^^ fe )||<A>;.A;.) 



it i' «1 

X < Ij I'. K 2 > < 

K12 K'12 fc 



S S (Tl 

S S (7 2 

t 0"12 cri 2 k 



> . 



(53) 



(4) Distributions a(3f3a, (3aa(3, j(3a-f, fl^ja, a-fSP, ja(35, (35ja, 5/3a^, a5jf3, 5a(3^, (3"/5a, -ff3a5 (cases 
4, 5, 13, 14, 28, 30, 33, 34, 37, 38, 41, 42 from Table 8) (Figure 1, A 4 , A s ): 



A 4 = 



E 



9 (niXi,rijXj,n' i X' i ,n' j X' j ,'5) A 8 , 
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Kl2d"l2^ 12 (J i2 
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eimXunjX^n^n'jX'^E) = | ^^+^+1',+^+^+^ x 
x [«i2,a 12 ,4 2 ,ai 2 ] 1/2 (n^XAlg^^'^^ln'MXr) x 



U l'i «i 

x < K 2 > < 

K12 K'12 fc 



S S (Tl 

S S cr 2 > . 
k CT12 cr' 12 fe ^ 



(55) 



(5) Distributions (3aaa, aflota (cases 7, 8 from Table 8) (Figure 1, A 5 ): 

A 5 = ^2 (n i X i ,n j X j ,ri i X' i ,n' j X' j ,'E.) > 



a (; aS ) x [ 5 (J„-) x 5 (J„-)] 



-1 (fefe) 



(56) 
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where 



Q(n i Xi,n j X j ,n' i X' i ,n' :j X' j ,'S) =6' (n i X i ,n j X j ,n' i X' i ,n' j X' j ,'E) = 
= (n a X a ,npXp,'E) . 
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and when G(T) = G(a(3aa), then 
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(6) Distributions PPfia, j3f3af3 (cases 9, 10 from Table 8) (Figure 1, A 5 ): 

A 5 = ^2 &(n i Xi,n j Xj,n' i X' i ,n' j X' j ,'Z) x 

Kl2CTl2^i2 (T 12 



[a( l f )S '> x o ( '" s) 



(K12CT12) 



X a('° s ) 



(fefe) 



where 



(niAi,njAj,n<A<,n^,S) = 6' (nA,nA, n i A i> ra j A j> s ) 
= (n a X a ,npXp,E) . 



When G(T) = G((3f3(3a): 
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and when G(T) = G{/3(3a(3): 



x [Ki2,a 12 ] 1/2 {np\pnp\ y^ k ^ k )\\n a \ a np\ p ) > 

l a Ifj k' 12 || s s &12 



X < 



Kl K 2 k > < C7l £72 }• X 

i/3 '/3 K 12 J ( S S (Ti2 



I K12 fc ^ J I 0-12 s J 

The final analytical expressions for diagram Ag appearing in (50), diagram Aj in (52) and diagram As in 
(54), are obtained after the following graphical transformations: 

(i) The second quantization operators are interchanged, until (from left to right) first come the operators 
acting upon shell a , then correspondingly upon (3, 7, <5. 

(ii) The generalized Clebsch-Gordan coefficient is transformed to match the order of operators. This is 
performed by changing the order of angular momenta coupling at some of the nodes 1, 2, 3 (Figure 1, 
A 6 , A7 and A%). 

We write down immediately the algebraic expressions for diagrams A e , A 7 and As (of Figure 1) after 
transformation (i) and (ii) by applying usual generalized graphical technique (see Gaigalas et al 1985). Also 
we have to notice that 6' is equal to with a phase factor, which is found by transforming the diagram of 
the tensorial structure according to the rules (i) and (ii). Only in cases 1 and 7-10 (see Table 8) is 9' = 9, 
because there is no need to transform the tensorial structure. 



As an example, let us consider in particular the case where operator aj acts upon the first shell niAi, 
operator ai acts upon the second shell n 2 A 2 , and operators at,, aj, act upon the third shell n 3 A 3 (see Eq. 
(1)). This is the distribution 18 in Tabic 8. We obtain the algebraic expression for distribution fiajj from 
(52). The two-particle operator for this distribution can be represented by diagram B\ which is proportional 
to its tensorial part (diagram B 2 ) as (Figure 2 B\, B 2 ): 



B x = ^ <d (n 2 \ 2 ,n 1 \ 1 ,n 3 \ 3 ,n 3 \ 3 ,~,) B 2 , (64) 

where 

9 (n 2 A 2 , mAi, n 3 A 3 , n 3 A 3 , S) = \ (-l) fe ~ p+1 [k 12 , ai 2l k[ 2 , a' 12 ] 1/2 x 
x (n 2 A 2 n 1 Ai||o( KlK2fe ' CTlff2fc )||n 3 A 3 a 3 A 3 ) x 

h h Ki 1 ( s s o\ 
x < h l 3 k 2 > < s s 02 > ■ (65) 
Kl2 k[ 2 k J [ CTi 2 cr' 12 k 

We use (i) for diagram B 2 (Figure 2, B 2 ) then. As in the expression (64) the order of second quantization 
operators is a^ 23 ^^ 1 ^^ 13 ^^ 13 ^, so we change it according to (i) and obtain (Figure 2 B\, B 3 ): 



B\ = - ^2 Q (n 2 X 2} n 1 X 1 ,n 3 X 3 ,n 3 X 3 ,^) B 3 . 

Ki2Cri2^' 1 2 (J i2 



(66) 
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Figure 2: Diagrams for distribution a^ 2 '^ a^ llS ^a,^ l3S ^cS l3S \ Diagram B\ represents a two-particle operator 
G. Diagrams _B 2 , B3 and B 4 represent graphical transformations. Diagram _B 2 represents the tensorial 
part of the two-particle operator B\ before transformations, diagram B3 represents this tensorial part after 
transformation i), and diagram B4 represents it after transformation ii). 
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We use (ii) for diagram B 3 then. We change the sign at the node 1 to obtain finally (Figure 2 B\, B4): 

'y ' ^_^h+l2+2s— K12— CT12+I x 

Kl2^12^i2 (J 12 

x6 (n 2 A 2 , niAi, n 3 A 3 , n 3 A 3 , S) B 4 = 
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(67) 



where 
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All the graphical transformations are made and the correspondence between angular momentum diagrams 
B 2 , B 3 and B4 in Figure 2 and algebraic expressions is defined according to the graphical approach of Gaigalas 
et al (1985). 

Now using (67) and (68) we can write down the irreducible tensorial form for the Coulomb operator with 
tensorial structure K\ = k 2 =k, a\ = a 2 = 0, k = and the two-electron submatrix element: 



(n 2 \ 2 n 1 \ 1 \\g { Couiomb I |n 3 A 3 n 3 A 3 ) = 
= 2 [fc] 1/2 (Z 2 ||C( fe )||Z 3 ) {h\\CW\\l 3 ) R k (n 2 l 2 n 3 l 3 , ni hn 3 l 3 ) . 



(69) 



From (68) we have: 



^'coulomb (n 2 A 2 ,mAi,n 3 A 3 ,n 3 A 3 , S) 

= \ E (-1)' 

x2[fc] 1/2 (Z 2 ||C( fc )||Z 3 ) (h\\C^\\l 3 ) R k (n 2 W 3 ,niZm 3 / 3 )x 
?2 '3 k 
h h k 
k 12 k' 12 



[ki 2 ,o-i 2 ,k 12 ,ct 12 J x 



s s 
s s 



cti 2 a' 12 



2+i;?-o-i2+fc r„ „ il/2 

[Kl2,Cri 2 J X 



^12<7 12 

x 2 ||C( fe )||Z 3 ) (Z 1 ||C( fe )||Z 3 )i? fe (nafans/s.niZmais) 



?2 ^3 
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(70) 



From (67), by (70), we finally obtain the following expression, in second quantized form, for the Coulomb 
operator for the case a = 1, /3 = 2 and 7 = 3: 
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(71) 



This kind of operator needs to be calculated when we are considering, for example, the matrix element 
(3s 2 3d 2 L 1 S 1 Q 1 L 2 S 2 Q 2 LS\ \H Coulomb \ \3s3d3p 2 L[S[Q[L' 2 S 2 Q' 2 L[ 2 S' 12 L' 3 S 3 Q' 3 L'S') . 
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